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Abstract 

Fermion iV-loops with an arbitrary number of density vertices N > d + 1 in d 
spatial dimensions can be expressed as a linear combination of {d + l)-loops with 
coefficients that are rational functions of external momentum and energy variables. 
A theorem on symmetrized products then implies that divergencies of single loops for 
low energy and small momenta cancel each other when loops with permuted external 
variables are summed. We apply these results to the one-dimensional Fermi gas, 
where an explicit formula for arbitrary iVdoops can be derived. The symmetrized 
iV-loop, which describes the dynamical iV-point density correlations of the ID Fermi 
gas, does not diverge for low energies and small momenta. We derive the precise 
scaling behavior of the symmetrized ./Vdoop in various important infrared limits. 

KEYWORDS: Fermi systems, Feynman amplitudes, density correlations, surface 
fluctuations 



1. Introduction 

The properties of fermion loops with density vertices (see Fig. 1) play a role in the the- 
ory of Fermi systems and various other problems in statistical mechanics. Symmetrized 
loops, obtained by summing all permutations of the N external energy-momentum vari- 
ables of a single iV-loop, describe dynamical iV-point density correlations of a (non- 
interacting) Fermi gas. Single loops have no direct physical meaning (for N > 2), but 
contribute as subdiagrams of Feynman diagrams in the perturbation expansion of inter- 
acting Fermi systems. Symmetrized loops appear as integral kernels in effective actions 
for interacting Fermi systems, where fermionic degrees of freedom have been eliminated 
in favor of collective density fluctuations JO. The behavior of symmetrized loops for small 
energy and momentum variables is particularly important for Fermi systems with long- 
range interactions, whose Fourier transform is singular for small energy and momentum 
transfers |2|, ||. 

Besides their relevance for interacting electron systems and other fermionic systems 
in nature, the theory of Fermi systems has also a bearing on various problems in clas- 
sical statistical mechanics, which can be mapped to an effective Fermi system (gas or 
interacting). For example, the statistical mechanics of directed lines in two dimensions 
can be mapped to the quantum mechanics of fermions in one spatial dimension Q. This 
mapping has been exploited extensively to study fluctuations of crystal surfaces |5|, [| . 

The 2-loop, corresponding to the 2-point density correlation function has been com- 
puted long ago in one, two, and three dimensions 0. Recently, Feldman et al. || have 
obtained an exact expression for the iV-loop with arbitrary energy and momentum vari- 
ables in two dimensions. We have evaluated that expression explicitly and analyzed the 
small energy-momentum limit of the symmetrized loops, showing in particular that in- 
frared divergencies of single loops cancel completely in the sum over permutations |§. 

Most recently, Wagner pll has published a reduction formula for fermion loops in the 
static case, where all energy variables are set zero. This formula reduces the iV-loop for 
a (i-dimensional Fermi system to a linear combination of (d + l)-loops, with coefficients 
that are rational functions of the momenta. In this work we point out that Wagner's 
formula and derivation can be easily extended to the case of finite energy variables (Sec. 
3). In the two-dimensional case, the possibility of such an extension is evident from the 
exact expression for iV-loops ||. The small energy-momentum behavior of symmetrized 
iV-loops can be analyzed by applying a theorem on symmetrized products derived in our 
work on two-dimensional systems [f|, which we formulate for the general <i-dimensional 
case in Sec. 4. We apply the reduction formula to a one-dimensional system, where the 
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iV-loop can be expressed in terms the 2-loop, which is very easy to compute (Sec. 5). We 
finally compute the infrared scaling behavior of symmetrized iV-loops in a one-dimensional 
Fermi system. 



2. Loops 

The amplitude of the iV-loop with density vertices, represented by the Feynman dia- 
gram in Fig. 1, is given by 

f d d k ( dk N 

U N ( qi ,...,q N ) = I N ( Pl ,...,p N ) = jj^y J^- TlGoik-pj) (1) 

at temperature zero. Here k = (k ,k), qj = (qj ,qj), and pj = (pj ,pj) are (d + 1)- 
dimensional energy- momentum vectors. We use natural units, i.e. h = 1. The variables 
qj and pj are related by the linear transformation 

Qj=Pj+i-Pj> 3 = 1,'-,N (2) 

where Pn+i = Pi- Energy and momentum conservation at all vertices yields the restriction 
qi + . . . + qN = 0. The variables qi, . . . , q^ fix pi, . . . ,p N only up to a constant shift 
Pj i — ► pj + p. Setting pi — 0, one gets 

P2 = qi 

P3 = qi + Q2 

Pn = Q1 + Q2 + ■■■ + Qn-i ■ (3) 
We use the imaginary time representation, with a non-interacting propagator 

G (k) = i (4) 

ik - (e k - /i) 

where e k is the dispersion relation and /i the chemical potential of the system. For a 
continuum (not lattice) Fermi system the dispersion relation is e k = k 2 /2m, where m is 
the fermion mass. The fco-integral in Eq. ([!]) can be easily carried out using the residue 
theorem; one obtains HI 
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Figure 1: The N-loop with its energy-momentum labels. 



where /y(k) = i(p i0 - p jQ ) + e k _ Pi - e k _ Pj .. 

The 2-loop n 2 (g, —q) = 14(g) is known as polarization insertion or particle- hole bubble, 
and has a direct physical meaning: 11(g) is the dynamical density-density correlation 
function of a non-interacting Fermi system For N > 2, the iV-loop is not a physical 
quantity, but the symmetrized iV-loop 

1 



n^Oi, ...,qN)=S H N (qi, • • • , gjv; 



An 



J2U N (q P1 , . . .,q PN ) 



(6) 



where the symmetrization operator S imposes summation over all permutations of q±, . . . 
qN, is proportional to the (connected) dynamical iV-point density correlation function: 



(p(qi),...,p(q N ))c 



ir-\N-l)\U%( qi ,...,q N ) 



(7) 



Here p(q) is the Fourier transform of the particle density operator. Eq. (7) is easily 
verified by applying Wick's theorem ||11|| . Note that Wick's theorem yields a sum of 
(iV — 1)! distinct loops with non-equivalent permutations of q±, . . . , q^, while the sum in 
Eq. (6) includes cyclic permutations which produce iV equivalent copies of each loop. 



3. Reduction formula 

We now state the reduction formula that reduces the iV-loop for a <i-dimensional 
system with iV > d + 1 to a linear combination of (d + l)-loops with coefficients that 
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are explicitly computable rational functions of momentum and energy variables. This 



formula is a straightforward generalization of a result derived recently by Wagner |10[ for 
the static case pjo = 0. 

Let pi, . . . ,p N be such that for each tupel of integers j = (ji, . . . , jd+i) with 1 < jx < 
■ ■ ■ < jd+x — N, the complex d-dimensional vectors d-> determined by the linear equations 

fnjM*) = i(p jl0 ~ P Jr0 ) + ^ (Pi - P 2 jr ) + ^ (P> - Ph) • d j = (8) 



for r = 2, . . . , d + 1 are well-defined and unique. Suppose that for n — 1, 
n 3x, ■ ■ ■ ,jd+x the numbers 



N with 



fl ■= /jVn(d J ) = i(p jr0 - Pno) + ^ (P 2 jr 



Pn) + — iPn ~ Pjr] 



are non-zero. Then 



(9) 



I N (px, ...,Pn) 



E 



i<n<...<j d+l <N 



N 

n 

n=l 



1 



I d+1 (p h , . . . ,Pj ) 



(10) 



Note that the numbers / Jrn (dj) with r 



d + 1 are all equal, as a consequence of 



Eq. (8). The vector d-> is uniquely defined if the vectors p Jr — p^, where r = 2, . . . , d + 
1, are linearly independent. The real part of d^ is the center of the uniquely defined 
circumscribing sphere through the points p JX , . . . , Pj d+1 in <i-dimensional euclidean space. 
In contrast to d-^, the numbers fl are invariant under a shift pj i— > pj +p and can thus be 
expressed in terms of the variables qx, . . . , qN- 

The proof of the above reduction formula, a simple generalization of the proof given 
by Wagner [TO] for the static case, is presented in the Appendix. 



4. Symmetrized products 

Symmetrized loops are obtained by summing over all permutations of external energy- 
momentum variables q±,...,qN as in Eq. (6). Up to a trivial constant, symmetrized 
iV-loops are the connected iV-point density correlation functions of the Fermi gas. The 
behavior of these functions in the infrared limit qj — > determines the long-distance (in 
space and time) density correlations, and is a crucial ingredient for power-counting of 
contributions to effective actions for collective density fluctuations. We will consider two 
important scaling limits: 

i) small energy-momentum limit liniA^o n^(Agi, . . . , XqN) , 
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ii) dynamical limit lim A ^ n^[(gi , Aqi), . . . , {q N o, Aqjv)] • 

Single AMoops diverge generally (for almost all choices of q\, . . . , q^) as X 2 ~ N in the small 
energy- momentum limit, which is what one would expect from simple power-counting 
applied to the integral (1). A notable exception is the socalled static limit, where the 
momenta tend to zero after all energy variables qjo have vanished. In that case one 
obtains a unique finite limit ITv — > ^N-iy. de N ~ 2 L =Ai ? where D(e) is the density of 
states [1^]. In the following we will show that systematic cancellations occur in the sum 



over permutations in the general small energy-momentum limit and also in the dynamical 
limit. 

The factor multiplying the (d + l)-loops in the reduction formula can be written as 

N I d+1 

**■•= n Tr = n^ (ii) 

n=l Jn r =l 

"#ii,...j'd+i 



where 



qj r+1 -l) = Un=j r +1 ji f O r jr+l > Jr + l 



r 




for j r+1 = j r + l 

Here jd+2 = ji, i-e. for r = d + 1 the index n runs from jd+\ + 1 to iV and then from 1 
to ji — 1. Note that depends also on differences of the energy-momentum variables 
Pj x , . . . ,Pj d+1 , besides the explicitly written arguments. As a product of M r = j r+1 —j r — l 
factors {fl)~ l , F* diverges as X~ Mr in the small energy-momentum limit, since each 
vanishes linearly. We define a symmetrized product 

S r (k h ..., k Mr +i) = ^ M - ^, H(kpi, • • • , k P{Mr+1) ) (13) 

where all permutations of k\, . . . , ku r +\ are summed. According to the following theorem, 
the symmetrized product S\ can be expressed such that the cancellations of singularities 
in the infrared limit become obvious. 

Factorization theorem: The symmetrized product S\ can be written as n^pni times a 
sum over fractions with numerators 

(k CT1 • k a[ ) (k CT2 • K> ) {K Mr ■ K> Ur ) (14) 

where cij ^ a\ and M r = j r+ % —j r — l, and products of 2M r functions f s as denominators. 
The functions / J have the form 

f(p,p) = KPo - p' ) + ^(P 2 - P 2 ) + (P - P) " d j (15) 
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where p = pj r and p' = pj r + (partial sum of hi, ... , kM r +i)- 111 each numerator, each 
momentum variable k\, . . . , &M r +i appears at least once as a factor in one of the scalar 
products. 

For example, in the simplest case M r = 1 one obtains 

k 2 ) + F}(k 2 , h) = — ™l k ^y -TTT (16) 

P \Pjr , Pjr [P jr , P jr + fo) 

The factorization theorem has been derived recently |J in the context of two-dimensional 
systems. The proof provides a concrete algorithm leading to the factorized expression. 
Since the algorithm is actually independent of the dimensionality of the system, we will 
not repeat the derivation here. 

The infrared scaling behavior of S-j, follows directly: 

i) Sj. is finite (of order one) and real in the small energy- momentum limit. 

ii) Si vanishes as \ 2Mr in the dynamical limit. 

To see this, note that the functions ft (p, p') vanish linearly in the small energy-momentum 
limit, and are purely imaginary to leading order in A, while they remain finite in the 
dynamical limit. 

The symmetrized product is thus much smaller for small energy and momentum vari- 
ables than each single term, namely by a factor X Mr in the small energy- momentum limit, 
and even by a factor X 2Mr in the dynamical limit. This result holds in any dimension d. 



5. One-dimensional systems 

We now apply the general results from Sees. 3 and 4 to one- dimensional systems |13|] , 
where particularly simple expressions can be obtained. We consider first single, then 
symmetrized loops. 

A) Single Loops 

In one dimension, the reduction formula (10) reduces TV-loops to linear combinations 
of 2-loops: 



I N (p h ...,p N ) = 
where d? is given explicitly by 



31,32 
\<31<32< N 



N 



n 7 

n=i Jn 

,32 



(17) 



2 Pjii -Phi 



and 

fn = ~W—(jPnl - Phl)(Pnl ~ Phi) + «(PjiO ~ Pno) + i(Pnl ~ Phi) — _ — (19) 

Here p n i and p Jr i are the one-dimensional momentum components of the energy-momen- 
tum vectors p n = (p n o,Pni) and pj r = (pj r o,Pj r i), respectively. The 2-loop can be com- 
puted very easily, the result being 



HPjnPh 



in 



T Phi ~ Phi 



log 



kp — acji 



J2 



kp + a 



where 



1 / \ , Pho ~ Pho 

a hh = 9 feii - Phi) + %m - — — 

z Piii Pj 2 i 



(20) 



(21) 



We have thus obtained an explicit expression in terms of elementary functions for iV-loops 
in one dimension. One may easily perform an analytic continuation to real (instead of 
imaginary) energy variables, ipjo i— » ej, in the above expressions to analyze, for example, 
the non-linear dynamical density response of the Fermi gas. 

In the zero energy limit pjo — > one obtains the simple result 



lim I N {pi, 

i=i,...,N 



E 



01,32 
l<j 1 <j 2 <N 



,PNj 



N 

S (Pnl-Phl)(Pnl-Phl) 

™#J1 ,32 



—2m 



rn 



n(Phi-Ph 



log 



2k F - (Phi-Phi) 



2k F + (Phi-Ph 



(22) 



Note that the above expression has a finite limit for pji — > 0, although each contribution 
to the sum diverges. 



B) Symmetrized Loops 

It is well known that for a linearized dispersion relation = f^G^I — kp), as in 
the one-dimensional Luttinger model, the symmetrized iV-loop U.fr(qi, ■ ■ ■ , <7at) vanishes 



identically for N > 2 even for finite qj with sufficiently small momenta qji ||14|| . We now 
analyze the infrared behavior of symmetrized iV-loops in a one- dimensional system with 
the usual quadratic dispersion relation. Symmetrizing the reduction formula, we can write 
symmetrized loops as 



gjv) = S Yl s i s 2 J 2 (Pii > Pi2 



(23) 



31,32 
1<J1<32< N 
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where S is the symmetrization operator introduced in Sec. 2 and S{ and S2 are the 
symmetrized products defined in Sec. 4. Note that first symmetrizing partially (with 
respect to a subset of variables, as in the products S-j.) and then completely (by applying 
S) yields the same result as symmetrizing everything just once. 

We can now easily derive the scaling behavior of IT^ in the small energy-momentum 
and dynamical limit, respectively. The 2-loop U(qi) = n 2 (gi,— q\) = hi^^qi) tends to 
the finite value 

n ( A ^i) - 1 _i_ r ), ™ ( 24 ) 

in the small energy-momentum limit and vanishes quadratically as 

n(g 10 ,A gil )--— ^ A 2 (25) 

in the dynamical limit, where v F = k F /m is the Fermi velocity. The same behavior is 
found for the 2-loop with a linearized e&. Since S{ and 5*2 are both finite in the small 
energy momentum limit, the symmetrized iV-loop remains finite, too: 

n s N (\ qi ,...,\q N ) = 0(l) for A-0. (26) 

Only in the static case qjo = each single loop Iljv has a finite limit for qji — > 0, while in 
general the above result is due to systematic cancellations of infrared divergencies. In the 
dynamical limit the product SfS^ vanishes as \ 2M i+ 2M 2 where M\ + M2 = N — 2, such 
that 

U s N [(q w ,Xqii),...,(qN ,Xqm)]=O(X 2N - 2 ) for A^O. (27) 
The same scaling behavior has been found previously for two-dimensional systems O. 



6. Conclusion 

We have derived a formula that reduces the evaluation of fermion loops with N density 
vertices in d dimensions to the computation of loops with only d + 1 vertices. This was 
obtained by a straightforward extension of a recent result by Wagner for the zero 
energy limit to arbitrary energy variables. Using a theorem about symmetrized products, 
we have shown that infrared divergencies of single loops cancel to a large extent when 
permutations of external energy-momentum variables are summed. The symmetrized 
iV-loop, which is proportional to the iV-point density correlation function of the Fermi 
gas, is thus generally much smaller in the infrared limit than unsymmetrized loops. For 
one-dimensional systems, we have obtained an explicit expression for arbitrary iV-loops 
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in terms of elementary functions of the energy-momentum variables. We have shown that 
symmetrized loops do not diverge for low energies and small momenta. In the dynamical 
limit, where momenta scale to zero at fixed energy variables, the symmetrized A-loop 
vanishes as the (2N — 2)th power of the scale parameter. 

We finally outline some applications of our results. 

Evaluation of Feynman diagrams: Analytical results for loops are of course useful for 
computing Feynman diagrams containing fermion loops as subdiagrams. The number of 
energy-momentum variables that remain to be integrated (analytically or numerically) 
is thus reduced. In particular, the mutual cancellation of contributions associated with 
different permutations of energy-momentum transfers entering a loop can be treated an- 
alytically, avoiding numerical "minus-sign" problems. 

Effective actions: Effective actions for interacting Fermi systems, where the fermionic 
degrees of freedom have been eliminated in favor of collective density fluctuations, contain 
symmetrized iV-loops as kernels [|IJ. A good control of the infrared behavior of these 
kernels is essential for assessing the relevance of non-Gaussian terms in the effective action, 
especially in the presence of long-range interactions. In one- dimensional systems one 
can use our results to compute the scaling dimensions of corrections to the leading low- 



energy behavior of Luttinger liquids []Toj| by analyzing the non-quadratic corrections in 
the bosonized action. 

Surface fluctuations: Some models of surface fluctuations lead to the statistical me- 
chanics of directed lines in two dimensions, which can be mapped to the quantum me- 
chanics of fermions in one spatial dimension || [|. Most recently, Prahofer and Spohn 
fl6|| have shown that the probability distribution of height fluctuations in such models 
is Gaussian at long distances on the surface. For this result it was enough to establish 
that symmetrized iV-loops in the associated Fermi system are less singular than the naive 
power-counting estimate. Our result Eq. (26) yields the precise scaling dimension of non- 
Gaussian terms, and implies in particular that high order corrections vanish very rapidly 
at long distances. 
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Appendix A: Proof of reduction formula 

Following Wagner's [HJ derivation for the static case, we prove the reduction formula 
(10) by applying the following many-dimensional version of Lagrange's interpolation for- 
mula: 

Lemma: Suppose that 1 < d + 1 < N and the (d + l)-dimensional complex vectors 
a 1; . . . , sl n are such that a^, . . . , a Jd+1 as well as (a^ — a n ), . . . , (a Jd+1 — a n ) are linearly in- 
dependent for pairwise different indices jx, . . . , jd+i,n G {1, . . . , N}. For j — (ji, ■ ■ ■ , jd+i) 
with 1 < ji < . . . < jd+i < N determine the complex (d + l)-dimensional vector z J by the 
system of linear equations a Jr -zj = 1 for r — 1, . . . , d+1. Then each complex homogeneous 
polynomial P(zq, z) of degree N — (d + 1) in the d + 2 variables zq, z = (z\, . . . , Zd+i) can 
be written as 



P(z ,z) 



E 



Pfl,z J ^ 



N 

n 



(2 -a n -z)det(a il ,...,a id+1 ) 



l<31<-<i'd+l<JV 





( 1 


1 . 




det 










\ a„ 


a ii • 


• • a id+i / 



(28) 



where the vectors a 1; . . . , a^v enter the determinants as column vectors. 
For a proof, see Ref. fl7 |. 

We apply the above lemma to the polynomial P(zq, z) = and 

-i(ko-Pno) + £ P71 
(k-2p n )/^2^ 

where £ p = p 2 / (2m) — /i. Since P(l, z^) = (z\) N ~( d+1 ^ and a Jr - z-> = 1 for r = 1, 
Cramer's rule yields 



(29) 
d + 1, 



TV 



N 



P(l,z j ) n det(a il ,...,a jd+1 ) = [J zj det(a ii; . . . , a jd+i ; 



n = l 
i#Jl,.--,Jd+l 



n = l 

n^Jl,---,j'd+l 



n det ( k-2 PJ1 k- 2p 

V2m 



n = l 



2i;i 



AT 

n 



n = l 



x d 

2/m J det (p ia - p^ , . . . , p id+1 - p j: ) 



(30) 



In the last step we have subtracted the first column of the determinant from all the others 
and then applied Laplace's theorem. We now evaluate the denominator in (|28|), 



D = det 



1 1 ... 1 

-i(ko-Pno) +t Pn -i(k -Pho) + £ Pjl ••• -i( k o~Pj d+1 o) + £ Pjd+1 
V (k- Pn )/V2^ (k-p n )/V2^i ... (k-p id+ J/V2m" J 



(31) 
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Subtracting the first column from all the others and applying Laplace's theorem yields 
D dot [ &" (0) J; d+A ° ] 1 (32) 

(Pn-Pji) ••• Vm (P« ~ Pid+1, 



Adding • (p n — pj r )/m to the r-th matrix element in the first row (adding thus multiples 
of the other rows to the first one) one obtains 



D = det 



„ (Pn — Pji) • • • Vm (P™ — Pid+i> 



2/ro) #det( 1 1 (33) 

' V P«"Pii ••• P«-Pi d+ i / 

Subtracting the first column from all others and applying Laplace's theorem once again 
one obtains 



D = {-y 2 / m ) fn det(p J2 -p ji; . . . , p jd+1 -p h ) (34) 

Eq. (bud and Ec i- (0) y ield 





det( 


, a ii 5 • • 


• ' a id+i) 




( 1 


1 




det | 














• • • a id+i / 



1 



n / /T" ; Jd+i N = n ^ (35) 

n=l I 1 1 1 \ n=l Jn 

n^il.— Jd+l det J'd+1 



We now set z = and z = (1, k/ i/2m), such that 

z - a n - z = z(fc„ - Pno) - ^k- Pn = Go 1 ^ - Pn) ( 36 ) 
With this choice of variables the above lemma thus yields the algebraic identity 

N 1 

1= E II -jjGffr-Pn) (37) 

l<il<-<J«J+l<iV «#il.-.i'd+l 

Multiplying this equation with Iljli Co(A; —pj) and integrating over one finally obtains 
the reduction formula Eq. (|T0|). 
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